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1 RESEARCH OBJECTIVES AND POTENTIAL IMPACT ON PROPULSION 

/ 

The rotating machines used in space engines; e.g., orbital transfer v e hic ' es - 
space shuttle main engines (SSME), operate at much higher speeds compared to 
those used in ground based engines or aircraft engines. The SSME oxygen 
turbopumps operate at about 30,000 rpm; and rotational speeds of space engines can 
be as hiah as 100 000 rpm. At these high speeds, conventional bearings, e.g., ball 
beahngs 9 are not reliable^ Consequently the main objective of this research program 
is to develop a highly reliable magnetic bearing system for space engine applications. 

In rebent Tears a number of researchers have developed magnetic bearing 
systems However, their main focus has been the applications to relatively low speed 
Ingines. Currently, NASA Lewis Research Center is developing magnetic bearings 
for SSME turbopumps. The control algorithms which have been used are based on 
either the proportional-intergral-derivative control (PID) approach or the linear 
quadratic (LQ) state space approach. These approaches lead to an acceptable 
performance only when the system model is accurately known, which is seldom true in 
practice. For example, the rotor eccentricity, which is a major source of vbration at 
high speeds, cannot be predicted accurately. Furthermore, the dynamics of a rotor 
shaft, which must be treated as a flexible system to model the elastic ro or shaft is 
infinite dimensional in theory and the controller can only be developed on the basis o 
a finite number of modes. Therefore, the development of the control system is further 
complicated by the possibility of closed loop system instability because of residua or 
uncontrolled modes, the so called spillover problem. Consequently, novel control 
algorithms for magnetic bearings are being developed to be robust to inevitable 
parametric uncertainties, external disturbances, spillover phenomenon and noise^ 
Also, as pointed out earlier, magnetic bearings must exhibit good performance at a 
speed over 30,000 rpm. This implies that the sampling penod available for the design 
of a digital control system has to be of the order of 0.5 milli-seconds. Therefore, 
feedback coefficients and other required controller parameters have to be computed 
off-line so that the on-line computational burden in extremely small. 

The development of the robust and real-time control algonthms is based on the 
sliding mode control theory. In this method, a dynamic system is made to move along 
a manifold of sliding hyperplanes to the origin of the state space. The number of 
sliding hyperplanes equals that of actuators. The sliding mode controller has two 
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parts: linear state feedback and nonlinear terms. The nonlinear terms guarantee that 
the system would reach the intersection of all sliding hyperplanes and remain on it 
when bounds on the errors in the system parameters and external disturbances are 
known. The linear part of the control drives the system to the origin of state space. 
Another important feature is that the controller parameters can be computed off-line 
Consequently, on-line computational burden is small. 

2. CURRENT STATUS AND RESULTS 

2.1 CONTROL ALGORITHM FOR RIGID ROTOR 

First, the flexibility of the rotor shaft has been ignored and the basic 
understanding of the bearing control system has been obtained. The equations of 
motion for the generic model of a rotor shaft (Figure 1) have been obtained using 

Newton and Euler laws. For small displacements, the linearized equations have the 
following form: 


X + PX = BU + VS, + D 


(1) 


where XeR and represents the radial displacements at locations of magnetic 
bearings. Elements of 4 x 1 control input vector V correspond to directions of the 
elements of C. The matrix P represents the gyroscopic effects. The 2x1 vector Se is 

simply (sin flt cos flt]T where D is the angular speed of the shaft. The vector D 
contains transient disturbance, and it is assumed that upper bound on each element, 
Qmax. is known. The uncertainty in the model due to the lack of precise knowledge of 
the rotor eccentricity is contained in the 4x2 matrix V. it is assumed that 



( 2 ) 


where Vjj is the estimate of Vjjand Fjj is the maximum error, which is considered to be 
known. 

The sliding mode control law is as follows: 


where 


(BU); = -XjXj - sinQt - cosiit - Kj satf ^ | 

S, = X, + x,x, 

Kj = F,i|sinftt| + F£|cosQt| + dmax+'ni 

4>i = Kj / *.j 
Xj>0 


(3) 


Here, rjj is a small positive number. Note that Sj=0 defines the sliding line (Figure 2) 
and 4>j is the time varying thickness of the boundary layer, which has been introduced 
around the sliding line to eliminate the chattering behavior. 
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A computer program has been developed to simulate the digital implementation 
of the control law, eq. (3). A representative steady state response of the controlled 
rotor is shown in Figure 3. The effects of X\ and sampling period Ts on the 
performance of the control system are being examined. 

A novel algorithm has been developed to determine the current in the magnetic 
bearing so that the control force required by equation (3) will be exactly achieved in 
spite of the nonlinear relationship among the magnetic force, coil current and the air 
gap. Using this algorithm, the level of maximum current will be estimated for the 
parameters of space engines. 

2.2. FLEXIBLE ROTOR DYNAMIC MODEL 

A mathematical model of the flexible rotor system has been formulated. The 
rotor system is modeled as a fixed-free axisymmetric shaft with an unbalanced disk 
inertia, and supported by two electro-magnetic bearings, or four independent actuators 
(Figure 1). The equations of motion and boundary conditions shown below are 
derived by applying the Hamilton’s principle. 

The equations of motion are: 

mx. tt+Elx.7777+(mffX,ti-lpQy.zzt-ltx,zztt)5(z-d) 

=Fi35(z-a)+F575(z-b)+pen 2 cosflt 5(z-d) (4) 

my, tt+Ely.zzzz+(mcjy.tt+lp^ix.zzt-lty.zztt)5(z-d) 

=F245(z-a)+F685( z -b)+P e ^ s ' n £M 5(z-d) (5) 

The boundary conditions are: 

x(0,t) = x, 2 (0.t) = y(O.t) = y,z(0,t) = 0 (6) 

x.zzO-.t) = x,222(l-.t) = y.zz(Ut) = y.zzz(Ut) = 0 (7) 


Here, 


Cl = shaft rotational speed E = elastic modules m^ = disk mass 

I = shaft second moment of area m = shaft mass per unit length L = shaft length 

p = disk unbalanced mass e = disk eccentricity 8 = Delta function 

Ip = disk polar moment of inertia It = disk transverse moment of inertia 

FI 3 = bearing force F1-F3 F24 = bearing force F2-F4 

F57 = bearing force F5-F7 F68 = bearing force F6-F8 

(),t = partial derivative with respect to time 

(),z = partial derivative with respect to space coordinate 


Space engine rotor models developed by other researchers will be examined. 
The important results from these previous analyses will be used to adjust and tune the 
above model to reflect the space engine rotor dynamic characteristics. The final model 
will be used as the basis for the dynamic analysis and control algorithm synthesis of 
flexible rotors. 


2.3 EXPERIMENTAL DESIGN 
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In parallel to the analytical efforts, an experiment has been planned and 
designed to validate the rotor model, implement the control algorithm, and verify the 
theoretical predictions. A rotor fixture from Bently Nevada Company (Figure 4) has 
been specified and purchased. The fixture consists of a shaft with lumped disk inertia. 
The unbalanced force of the rotor can be adjusted through insertion of weights into the 
disks. The shaft speed is controlled by an AC motor. Measurements of the transverse 
vibration of the rotor are performed by using non-contact displacement probes. Two 
magnetic bearings are being designed and fabricated by Magnetic Bearing 
Corporation. The rotor fixture is presently being modified to incorporate the magnetic 
bearings. A PC-based micro-processor control system has been specified and is 
being set up as the main controller for the control algorithm implementation. 

3. PROPOSED WORK FOR COMING YEAR 

3.1 CONTROL ALGORITHM FOR FLEXIBLE ROTOR 

The number of vibratory modes of a flexible rotor is infinite in theory and 
extremely large in practice. Since it is not practical to control vibration in all the 
modes, the controller will be designed on the basis of a finite number of modes, which 
is termed as controlled modes. The remaining higher frequency modes are called 
'residual modes'. Since magnetic bearings will provide control forces in four 
independent directions, the effective numbers of actuators will be less than the number 
of controlled modes, N. In this case, sliding hyperplanes will have to be taken as S = 
EX where S and X are respectively 4 x 1 and 2N x 1 vectors, and E is a full matrix 
unlike the situation for the rigid rotor. Furthermore, all the elements of the state vector 
X cannot be directly measured. Consequently, these states have to be obtained using 
an observer. The controller design will essentially involve an appropriate choice of the 
matrix E and the observer gain matrix. Using the singular value robustness tests, 
these matrices will be chosen such that the closed loop system is asymptotically stable 
in the presence of residual modes also. Steady state analyses will be performed to 
determine the influence of matrix E on the magnitude of vibration level and the control 
force. The objective of this analysis will be to determine an optimal E for given 
characteristics of magnetic bearings. 

The applications of model reference, sliding mode adaptive control technique 
will also be investigated. In this approach, the sensor output can be made to behave 
like a response of a reference model having desired damping ratio and natural 
frequency. The rotor response at those locations where sensors are not mounted will 
be investigated. 

Rotor parameters of space engines will be used in designing these control laws 
vanous implementation requirements such as power, current etc. will be estimated. 

3.2 EXPERIMENTAL VALIDATION OF THE CONTROL SYSTEMS 


Two phases of hardware development will be carried out as follows: 
Phase I - Test Stand Set Up and System Characterization 
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The rotor fixture, the magnetic bearings, the sensors and the controller will be 
integrated and a shakedown test will be performed. A series of tests will first be 
carried out to determine the magnetic bearings' dynamic characteristics. The rotor will 
be excited at various positions with impulsive, step, and periodic inputs, in order to 
characterize the structural dynamics. The results of these efforts will be used to 
validate the analytical model and modify the control law. 

Phase II - Control Implementation and Validation 

In the second phase of the experimental study, sensors and actuators for control 
purposes will be applied to the test stand according to the recommendations from the 
analytical work. The rigid rotor control algorithm will first be tested. System 
parameters, such as sensor locations, shaft length, shaft speeds, and rotor 
unbalanced force, will be varied systematically to examine the performance and 
robustness of the controller. The results from this phase will provide a quantitative 
measure of the efficacy of the proposed control strategy. 
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Figure 3: Steady State Response of the Controlled Rigid Rotor 
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